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In this paper we present a comprehensive formaUsm for dilepton production from the collision 
of two polarized spin-i hadrons by identifying the general angular distribution of the cross section 
in combination with a complete set of structure functions. The various structure functions are 
computed in the parton model approximation where we mainly consider the case when the trans- 
verse momentum of the dilepton pair is much smaller than its invariant mass. In this kinematical 
region dilepton production can be described in terms of transverse momentum dependent parton 
^\ , distributions. 
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I. INTRODUCTION 



' During the past decades dilepton production in high-energy hadron hadron collisions (the so-called Droll- Yan 

CO . (DY) process HI, [H) has played an important role in order to pin down parton distributions (PDFs) of hadrons. 

While the main focus was on PDFs of the nucleon, also information on the partonic structure of the pion was 
' (— I ', already obtained through Drell-Yan measurements. The crucial tool required for the extraction of PDFs is the 
pi i' QCD-factorization theorem [1, which applies if the invariant mass of the dilepton pair is sufficiently large. 

Experimentally, the Drell-Yan process is quite challenging because of the relatively low counting rates. On the 
other hand, from the theoretical point of view it is the cleanest hard hadron hadron scattering process. The fact 
^ I that no hadron is detected in the final state simplifies the proof of factorization in comparison to hadron hadron 
collisions with hadronic final states. This important point is one of the main reasons for the continued interest in 
the Drell-Yan reaction. 

Currently, not less than six programs for future Drell-Yan measurements are pursued. These plans comprise 
dilepton production in nucleon nucleon collisions (at RHIC 0], J-PARC (KEKljl, i], IHEP (Protvino) [13], 
I and at the JINR (Dubna) [11]), in antiproton nucleon collisions (at FAIR (GSI) [I4]), as well as in pion nucleon 
■ collisions (at COMPASS (CERN) 13]). Past measurements exclusively considered the unpolarized cross section, 
O^l I but all future programs are also aiming at polarization measurements. Including polarization of the incoming 
. hadrons opens up a variety of new opportunities for studying the strong interaction in both the perturbative and 
' the nonperturbative regime. Here we only mention the access to the transversity distribution of the nucleon (l^ . 
00 ! [H, HB, E [H, H |2l|, and to transverse momentum dependent parton distributions (TMDs). The TMDs 
' not only depend on the longitudinal momentum of a parton inside a hadron but also on its (intrinsic) transverse 
momentum and, in general, describe the strength of various intriguing spin-spin or spin-orbit correlations of the 
parton- hadron system (see Refs. [l^jlH, [U[2^ for more information on TMDs). 

In order to analyze upcoming data from polarized Drell-Yan measurements it is necessary to have a general and 
^ ■ concise formalism at hand. The main motivation for writing the present paper is to provide such a framework. 
To this end we decompose the hadronic tensor of the polarized Drell-Yan process in terms of 48 basis tensors 
which are multiplied by structure functions. We limit ourselves to photon exchange and do not consider weak 
interaction effects. To ensure electromagnetic gauge invariance of the hadronic tensor we make use of a projector 
method proposed in Ref . [11] . On the basis of the hadronic tensor we then write down the general structure of 
the angular distribution of the Drell-Yan process. This step is most conveniently done in a dilepton rest frame 
like the Collins-Soper frame [13]. In addition to our model-independent results we also consider the process in 
the parton model approximation, where we distinguish between two cases: (1) cross section integrated upon the 
transverse momentum qx of the dilepton pair; (2) cross section kept differential in qx and qx <C q, where q is the 
invariant mass of the dilepton pair. While in the former case one ends up with ordinary forward PDFs, in the 
latter TMDs enter in the parton model description and in a full QCD treatment (28l. [29l. [soj . 

In addition to our model-independent treatment we also consider the process in the parton model approximation 
by concentrating on the situation when the cross section is kept differential in the transverse momentum qx of the 
dilepton pair. In this case TMDs enter the parton model description as well as a full QCD treatment [28i. l29l[30t . 
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FIG. 1: Amplitude for dilepton production in parton model approximation. Both diagrams have to be taken into account. 
The spectator systems Xa and Xt of the two hadrons are not detected. 

Part of the results presented here were already given in the literature [13, Isil . Is^ . [33^ , and we comment on other 
work during the course of the manuscript. However, to the best of our knowledge, a complete formalism for the 
polarized Drell-Yan process has not been worked out before. 

The manuscript is organized as follows. In Section II we fix part of our notation and give the general form of the 
cross section in the one-photon exchange approximation. Section III contains the decomposition of the hadronic 
tensor in terms of basis tensors and structure functions, while in Section IV some discussion on reference frames 
is given. In Section V we present the general angular distribution of the polarized Drell-Yan process which can 
be derived from the results of Section III in a straightforward manner. Section VI contains the results for the 
structure functions in the parton model approximation. We conclude in Section VII. 



II. CROSS SECTION IN ONE-PHOTON EXCHANGE APPROXIMATION 



To be now specific we consider the dilepton production 

with {Pa,Sa) and (Pf,,S'b) denoting the 4-momenta and the spin vectors of the incoming hadrons. One has 
— M^, Pa ■ Sa ~ 0, — ~l, and corresponding relations for the second hadron. Throughout this work the 
mass of the leptons in the final state is neglected. We will sum over the helicities A, A' of the leptons. 

At large invariant mass q of the dilepton pair the process (H)) can approximately be described in the Drell-Yan 
model [l, 0| , which corresponds to the parton model approximation. According to this approach a quark from 
hadron Ha and an antiquark from hadron Hb (and vice versa) annihilate into a timelike virtual photon which 
subsequently decays into a lepton pair (see Fig. [T]).^ This means the process proceeds according to 

Ha + Hb^-f*{q) + x + 1+ + X , (2) 

where the 4-momentum of the virtual photon is given hy q ~ I + I' Note that the meaning of ^ remains valid 
if higher order QCD corrections are taken into account. 

In the one-photon exchange approximation the (frame-independent) cross section of the Drell-Yan process is 
given by 



da 



em 



dHdH' Fq^ 

where 



ip.W^^'' , (3) 



F = AJ{Pa-PbY-MlMl (4) 



^ As already mentioned we do not consider weak interaction effects. 

^ In our notation the symbol q describes both the 4^momentum of the virtual photon as well as the invariant mass ^/(p of the dilepton 
pair. This should, however, not lead to any confusion. 
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represents the flux of the incoming hadrons. If hadron masses are neglected one can write F ~ 2s ~ 2{Pa + PbY ■ 
The fine structure constant is related to the elementary charge through aem — /Air. In Eq. ([3|) the quantity 
L**" denotes the spin-averaged leptonic tensor, 

L'^- = ^ {u{l,\)Yv{l' ,X')){u{U\)Yv{l\>!)y ^aUi'-^ + I'^I'^ - ^g^A , (5) 



A, A' 



while 



W^''{Pa,Sa;Pb,Sb;q) = I d'^xe^'^-'' {Pa,Sa;Pt,Sb\J^^{0)J^,M\Pa,Sa;Pb,Sb) (6) 

is the hadronic tensor, which is determined by the electromagnetic current operator J^„. 

The tensor W'^'^ a priori is unknown and contains the information on the hadron structure. It has to fulfill 
certain constraints due to electromagnetic gauge invariance, parity, and hermiticity. In this order the constraints 
read 

q^.W''''{Pa,Sa;Pb,Sb;q) = q,W''''{Pa,Sa;Pb,Sb;q) = 0, (7) 
W^^iPa^Sa-^Pb.Sbiq) = W^,{Pa,-Sa;Pb,-Sb;q), (8) 

W'^-'iPa^SalPb^Sbiq) = [W''^{Pa,Sa;Pb,Sb;q)]* , (9) 

where the notation — for a generic 4-vector v is used. In Section III, by imposing the relations 

the hadronic tensor is decomposed into a set of 48 basis tensors multiplied by scalar functions (structure 
functions). In doing so the conditions ^ and ^ considerably reduce the number of allowed basis tensors, while 
the hermiticity constraint ^ implies that the structure functions are real. Note that time-reversal does not 
impose any constraint on the hadronic tensor, because this operation converts the two-particle hadronic in-state 
into a two-particle out-state, and both states are not related. In Section VI the hadronic tensor is considered in 
the parton model approximation. 

The angular distribution of the Drell-Yan cross section is most conveniently be considered in a dilepton rest 
frame like the CoUins-Soper frame ^273 or the Gottfried- Jackson frame j34j. In any dilepton rest frame, one can 
rewrite Eq. ^ according to 

- """^ ■ L^.W" , (10) 



Sq dn 2 F g4 

where the solid angle specifies the orientation of the leptons. In Section IV we elaborate a bit more on reference 
frames with the main focus on the center-of-mass frame (cm- frame) and the CoUins-Soper frame (CS-frame). 



III. HADRONIC TENSOR 



The total hadronic tensor can be decomposed into the unpolarized, single polarized (for hadron Ha and hadron 
Hb), and double polarized tensor according to 

^rf^^ ^ ^r^^'^ + ^M-' + ^r^'^ + 1^^^^ . (11) 

In the following we merely have to consider the symmetric part of W^'^ because the spin-averaged leptonic tensor 
in ^ is symmetric under the exchange fi ^ v. 



A. Unpolarized case 

Since the unpolarized tensor depends on the 4- vectors q'^, P^, and P^^ one can immediately write down the 
tensor basis 

"■u.2 — q q 1 
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= Pa Pa , 




Ka 


= PbPb , 






= 9' -Pa ^ 


-1 Pa ^ 




= q^'Pb ^ 


l-'Pt, 


"■u,7 




1 ^a 



(12) 

The expressions in ((T^ constitute a complete list of basis tensors being in accordance with the parity constraint ([H]). 
Therefore one can write in a first step 

7 

4=1 

where the structure functions Vu.i depend on the invariants Pa ■ q, Pb ■ q, and q^. 

So far we have not yet used the gauge invariance constraint ([7]) which, in fact, implies that not all of the Vu^i 
are independent. Contracting the tensor in (jl3p with the 4- momentum of the virtual photon and imposing ([7]) 
one readily finds 

= Vu,i+q^K2 + Pa-qVu,5 + Pb-qVufi^ 

= Pa-qK3 + q^Vu,5 + Pb-qK7, 

- Pb-qVuA + q'Vus + Pa-qKr- (14) 

These three relations follow because in W^'^qi, the terms proportional to q^ , P^, and P^^ must vanish separately. 
Now one can use (jl4p to eliminate three structure functions and consequently ends up with a hadronic tensor 
given by just four independent structure functions that are multiplied by four independent basis tensors. The 
explicit form of the basis tensors depends of course on which of the structure functions are eliminated. 

Though this procedure of implementing gauge invariance in principle is straightforward it gets rather cumber- 
some for single and double polarization because in those cases considerably more structure functions and basis 
tensors are involved. Therefore we resort to an alternative and very elegant method proposed in Ref. [1^ which 
makes use of projection operators. We define^ 

pM- = g'^- - , (15) 

q 



and let this operator act on the basis tensors in (|T2|) according to 

P'phZPa"- (16) 

Because of the property 



the tensors in p6p vanish for i — 2, 5, 6, while the remaining four nonzero tensors are gauge invariant by 
construction. This means that one arrives at the following final form of the unpolarized hadronic tensor: 
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with the four structure functions V,, , , and the tensor basis 



^uA - 9 ^2 

''«,2 ~ ^a ^a ^ 
''«,3 ~ ^b ^b 1 



C4 = Pa'Pb+P^Pb- (19) 



^ Note that the projection operator is not unique l26l . One can also define an operator involving the hadron momentum Pa or Pf,. 
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In Eq. we make use of the vectors 

which vanish upon contraction with q. Needless to say that the tensor in (|18p is frame-independent. The 
contraction with the leptonic tensor can therefore be performed in any frame. In the context of the parton 
model calculation in Section VI, for instance, this contraction is carried out in the cm-frame. The specific form 
of the tensor (IT5)l is by no means unique. Other sets of basis tensors can be found in the literature (see, e.g., 
Refs. [13, H^, Sfil and also [131), and it is straightforward to write down relations between different sets. Here we 
refrain from doing so because it does not give much further insight and, in addition, is not needed for the main 
purpose of this paper. We have discussed the unpolarized case in some detail in order to outline the procedure 
which is used in the following two subsections that are dealing with hadron polarization. 



B. Single polarized case 

Now we proceed to the case when one of the hadrons in the initial state is polarized. We first consider 
polarization of the hadron Ha, and then just quote the result for the case when iJf, is polarized. In order to 
construct a tensor basis we now have also the spin vector Sf^ at our disposal — in addition to the 4- momenta g'^, 
P^, Pj^ . Imposing the parity constraint ([S]) one finds the following list of tensors which are symmetric under the 
exchange ii v: 



[Sa-q, Sa-Pb} (e^'^"^^ q" + e''^^"^' qn , 
[Sa-q, Sa-Pb} (e^«^"^^ p: + e^'iP-P' Pl^) , 

[Sa-q, Sa-Pb} (e^'^"^^ + e''"''"''^ Pb) ' 

gMSagPa ^ ^•^S.qP. q^ ^ 
e^^S.P.P, q. ^ ^.S.P.P, q^. ^ 

flSaqPa pf I pVSaqPa 

tiSaqPb pf I pvSaqPb 

^tiSaPaPb pf ^ ^vSaPaPb 
^USaqPa pV I ^uSaqPa pM 

£ ^b ' ^ ' 





' "-0,7 




' "-0,9 


"a, 10' 




"a, 12' 


"a, 13 




"0,14 




LA"' 
"■a,15 




"■0,16 




"a,17 




"a, 18 




"■0,19 




"'a,20 




"■a,21 




"■a, 22 




"a, 23 



a ' 



^fiSaqPb piJ _|_ ^vSaqPb pfJ- 

^flSaPaPb pl^ _|_ ^uSaPaPb 

^M9PoP. ^.^^g-'qPoP. (21) 

To shorten the notation we have used abbreviations like e^aqPaPb — e^^^p^S^q" P^Pj^ . Note that the hadron spin 
vector can only appear linearly. It turns out that not all of the tensors hZ in (|2ip are independent of each other. 
The identity 

ga/3 ^iivpa „ ^aupa _|_ ^i//3 ^fiapa _|_ ^p,vacT _|_ ^pvpa (22) 

allows one to eliminate several out of the 23 tensors. To be explicit one finds 10 linearly independent relations 
between the tensors in ([21]) which may be written in the form 



,23 ' 

-2 
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2h^^, = -Pa ■ Pt + Ml -Pa-q K:19 , 

= K^s - Ml Ka, + Pa-P, K':,, -P,-q Ka, , 



= Ksa -Pb-q +Pa-q ~ 9 ^a.W ^ 

_ p p , pi. , ^2 .A"' - P . n h^"" 



16 ' 



"a, 6 


— "a, 9 


-MlhZ, + Pa-PbhZ,- 


-Pb-q K'^w : 


"a, 6 


— "a, 12 


- Pb- qKao + Pa- qKai 




"a, 7 


^ "-a, 11 


-MlK:,, + p,-p,hZs 


-Pb-q Ka9 


"a, 7 


= -Pa- 


■PbhZo + Mlh^.-P,- 


q K%2 ■ 



(23) 

On the basis of the relations in we choose to chminate the tensors hll^^^, . . . , ft-^ 23- 

FoUowing Eq. (|16p the projection operator P^"^ is now apphed to the remaining tensors in order to implement 
electromagnetic gauge invariance. This procedure provides, in a straightforward manner, the final form of the 
hadronic tensor for the case of single hadron polarization. One finds 

8 

W^a^'' = EC^^".- (24) 

2 = 1 

with the eight structure functions Va^i, and the tensor basis 

.IIV .fj.l' _ SaqPaPb j r.t^'^ pfi fjl^ f}!^ pV Tjfi , pi' pM 1 

''a,n ■ ■ ■ ' '-a, 4 — = |.9 ^2~ ' ^ ^b ^b : ^b + ^a^b ( ^ 

C5' Ce ^ [Sa-q.S^-P,^ P: + e^"^"^^ P^) , 

= {Sa-q, Sa- Pb] (e''«^"^^ Pb^ + e'^"^"^^ ■ (25) 



''a, 7' a, 8 

Here we used the 4- vectors and as given in (PO)) . Note that the first four tensors in (|25p correspond to 
the four tensors in (|19p for the unpolarized case, multiplied by the structure sSaqPaPb_ jg worthwhile pointing 
out the following: we have chosen to first remove redundant tensors in (|2ip by means of the identity ([52]) and 
then implemented gauge invariance. If one reverses these two steps one can obtain the same final result for the 
hadronic tensor. 

If the hadron Hi, is polarized one can now write immediately 



w^r = E<»^M, (26) 



with the eight structure functions Vh^i, and the tensor basis 



,fLV ,fiV 

''6,1' • ■ • ' ''b,4 


^ ^StqPiPa 


'■b,5' ''b,6 


= ^Sb-q, 


''bj' ''b,8 


= ^Sb-q, 



ui/ q q tdp- pv pi' I pi' p^* 

9 I ^a^a 7 ^b ^b J -^a-^b ' ^b 



,Ai^i/ 

.fiqPhPa pf _j_ ^lyqPbPa 

^JjqPbPa _^ ^<^qP,Pa pP) _ (27) 

In Ref. (SlJ the case of single hadron polarization for the Drell-Yan process was already considered. In that 
paper, however, the focus of the model-independent part was on the angular distribution of the cross section in 
the CS-frame rather than on the general form of the hadronic tensor. We will discuss the angular distribution of 
the cross section in Section V. 
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C. Double polarized case 



Eventually, we consider the situation when both hadrons in the initial state are polarized. In that case the 
basis tensors depend linearly on both Sa and Sb- A full set of tensors respecting the parity constraint ([8]) and 
being symmetric under the exchange ^ ^ v reads 



"■ab.l' ■ 



"-06,8' ■ 



"■06, 22' ■ 



"ah, 29' ■ 



"■ab,36^ ■ 



"ah, 39' ■ 



' "ah, 7 



' '*ah,14 



''oh, 15' ■ ■ • ' '"ah, 21 



' "ah, 28 



' "ah, 35 



' "ah, 38 



' "ah. 41 



"-abAl' ■ • ■ ' "-abAA 



' "ah, 47 



"-ab^S 



Sa ■ St {9'^- , q^'' , Pi^P: , P^P^ , q^P: + q^'PH , g^'P^T + <f ' 
Pi^Pb+P^P^}. 

Sa-qSt-q {g^" , q^q'' , P^^ P^ , P^^ P^' , q^P^ + q" P^^ , g^P,^ + q"" P^' , 

Sa-qSfPa {g""" , q'^q'' , P^^P^ , P^P,"^ , g^P^ + q" Pi^ , g'^P,'' + g^P^ ' 
5h • 9 5„ • n {g^'' , q^'q" , Pa^P: , P^Ph"^ ' g^P: + g'^P^^ ' g^'Ph'' + g^'Pf 

Pa'Pb+P^Pt}^ 

Sa -PbSb- Pa {9''" , q^q" , Pa^P: , P^^ Pb' , q^ P^ + g'^Pa" , q"" P'b + g-'Pf 

PUPI^ + p^pC], 
Sa ■ q {5,^^ + su^ , s'^p: + S'^P^ , Sl'P,^ + S^P,^} , 

Sb ■ q {s^'g" + s:q^^ , s^^p: + s:p^ , Si: pi: + s^pt} , 
Sa ■ Pb {s^q" + su^ , sf^p: + s'ip:: , ^rn^^ + s'^p^] > 

•Jh ■ -Ta I "Jq g + '-'al I "-"a ^a + '-'a ^a I "-"a ^h + '-'a ^h J ' 
^a^b + "^a^b ■ 



h ' 



1/ 

h ' 



(28) 



Like in the case of single hadron polarization not all 48 tensors in (|28|) are independent of each other. An explicit 
relation between a certain subset of the /i^^j can be found by means of the determinant identity [38i | 



gfj.l' gfia gfj.i3 ^ 7 ^ /J ^ 

gOia gOfi gOy gOS 

gP. gPa gP-p ^/S^ ^ ^ 

g^" g'^" g'^^ g''^ c/T^ 

gSU gSa g&P gSy gSS 



0. 



(29) 



Equation (|29p immediately implies 



D 



a3^6 i^a^b + ^b) ^b ^b — , 



(30) 



which allows one to eliminate exactly one out of the tensors in (pS)) . For the sake of symmetry we choose to 
eliminate the tensor /iQ^4g. Equation pop implies a relation of the type 



{g2 [{Pa ■ Pbf - MlMl] -2Pa-PbPa-qPb-q + Ml {Pb ■ qf + Ml {Pa ■ g)2 } K'^^ 



48 



(31) 



where the r.h.s. of is a linear combination of terms in which most of the h'^^b i ~ 1, ■ ■ • ,47) enter. We 
refrain from writing down this (rather lengthy) formula explicitly as it is not needed for the following discussion. 
We also mention that the determinant identity does not lead to any further relation between the /i^j j. 
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To implement gauge invariance we now apply, according to Eq. (|16p . the projection operator P^'^ to the tensors 
in (|28p . This procedure provides, in a straightforward manner, the final form of the hadronic tensor for the case 
of polarization of both hadrons. One finds 

28 

W^aT^EC.KM, (32) 
with the 28 structure functions Vat.i, and the tensor basis 

.flU .fj.!/ rt c J „A"' p/J TDf pf^ fjv I pM 1 



j-fJ^l^ j-A"^ _ C . n ■ n } ni^^ ^ F"^ P" P'^ P" Pt^ P" ^ P" P'^ 

'■ab,5' ■ ■ ■ ' '-06,8 — '-'a 1 "^b </ S 5 " n ^ i ^b ^b ^ + ^b f' 



,Atjy _ r, . „ c p J ^M'' _ '^^'^ PA' pt^ pf^ pf PM p'' P'^ pl^ K 
''ab,9^ ■ ■ ■ ' ^ab, 12 — "-"a </ Jfc J^a <| i/ ^2 ' J ^b -^b ; ^b + ^a^b f ' 



j-A''' j-A"' _ . P, ^ nA"' _ pAi pi' pA* pi' pA' pi^ P^ Pt^ K 

Vb.lS' ■ ■ ■ ' ^ab,lQ — ^b r-b<^g , , ^f, , + ^ , 



iA"' 4-A''' c P Q P J ^A"' ^ ^ pAi p'' pA^ pi' PA" pf 1 pi' pM 

'■ab, 17' ■ ■ ■ ' ''ab, 20 ~ "Ja ' ^b ^b ' S </ 2 ' i ^b ^b 7 -^^fc + -^^a -f^b f ' 



,^1/ j./ii' _ Q „ J CAi pi' I CI' pAi QAi pi' I CI' pA' 
''afc,21' ^Qb,22 — '^a 1 S'^b 'T' "^b 1 '^b^b '^b ^b f' 

,^11/ ,/ji' _ . n J ^ Ai p'^ I CI' PAi CAi p'^ + P'^ , 
■-ab, 23' ^ab, 24 ~ y ^ "^a 1 '-'a^a ' '-'a ^b 1 "^a^b T' 

,^11/ ,/ji' Q p J cA* pi' I CI' PA* CA* pi' I Ci' pAi I 

''ab,25' ''ab,26 " ' ^b \Jb + -^b ' -^b ^b + -^b -^^b r ' 

,/^i/ ,/ji' _ Q, . P J QAi pi' I CI' PAi CAi pi' I Ci' pAi l /ooA 

''ab, 27' '-ab, 28 ~ ^a S ^a ^a +'-'a-'a ' '-'a ^b +'-'a-'b r ' V'J'Jy 

Here we used the 4- vectors and Pj^ as given in (|20p. The vectors S*^ and S*^ are defined accordingly, i.e., 

cAi _ cAi ' QQ^ Qfj. _ Qfj, Sb ■ qq^ , , 

qZ u u qz 

Note that the first 20 tensors in ([55]) correspond to the four tensors in (fTO)) for the unpolarized case, multiplied by 
certain scalar products containing the spin vectors of the hadrons. We again emphasize the crucial importance of 
the relation ([5n| . Without this identity the final form of the hadronic tensor would have 29 rather than 28 basis 
elements. 

To the best of our knowledge the general structure of the hadronic tensor for the double polarized Drell-Yan 
process is a new result. Though the double polarized case was already investigated in Ref. _T4,], this was only 
done for the specific cases qx — Q and cross section integrated upon qx- In those cases seven basis tensors can be 
identified. 

D. Identical hadrons 

If both hadrons in the initial state are identical — as is the case, e.g., for proton- proton DY — the total 
hadronic tensor in Eq. (fTTj) has to satisfy the symmetry relation 

W^'''{Pa, Sa, Pb, Sb-, q) = W^^'iPb, Sb-, Pa, Sa;q) ■ (35) 

This immediately implies that eight out of the 48 structure functions are symmetric when exchanging the momenta 
Pa and Pb, 

K,i(&, a) = K,i(a, 6) , Vu,A{b,a) = K,4(a, 6) , 

Vab,i(&, a) = Kb,i(a, 6) , Kb,4(&, a) = Kb,4(a, 6) , ^6,5(6, a) = 146,5(0, 6) , (36) 
Kb,8(&,a) = Kb,8(a, 6) , Vab,i7(&,a) = T4b,i7(a, 6) , V;b,2o(&,a) = Vab,2aia,b) , 
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where, e.g., the first relation in is a shorthand of 

VuAiPb -q.Pa- g, g') = K.iC/'a • <z, n • <z, q^) . (37) 

Because of the symmetry property it is sufficient to fcnow the structure functions in (j36p for just half of the allowed 
parameter space. The remaining 40 structure functions fulfil the relations 



Kt,3(&, a) — 


Ki,2(a, h) , 










*^ 0.1 ^) 




Vh o(b a) — 


V -lia b) 




V 'y(a b) 






VbAb,o,) = 


Va.,7{a,b) , 


Vbs{b,a) = 


Va,8{'^,b) , 


Vb.7{h,a) = 


Va,b{a,b) , 


V6,8(&,a) = 










= VabA°-^b) , 


Kfcj(&,a) = 


= Vabs{a,b) , 


Vabssib, a) 


= Vab,9ia,b) , 


VabM{b,a) 


= K6,ii(a,fc) , 


Vab,i5{b, a) 


= K6,io(a,fe) , 


Vabseib, a) 


= K6,i2(a,fe) , 


Vab,iSi{b,a) 




Vab,23{b, a) 


= Vab,22{a,b) , 


VabMib, a) 


= Kb,2i(a,^) , 


Vab,27{b,a) 


= V'ab,26(a,fo) , 


Vab,2&{b, a) 


= Kf,,25(a,^) • 







(38) 



For instance the first relation in psp implies that if one knows the structure function Vu^2 for the entire parameter 
space one also knows Vu,3- 



IV. REFERENCE FRAMES 



So far our treatment is frame-independent. If, however, one wants to write down the general form of the angular 
distribution of the cross section — as we are going to do in Section V — one has to specify the reference frame. 
Moreover, the parton model calculation of the hadronic tensor, carried out in Section VI, is naturally performed 
in the c?7i-frame. Therefore, in the following we will consider both the cm-frame and the CS-frame [27], which 
is a particular dilepton rest frame. In general, the angular distribution of the cross section is most conveniently 
given in a dilepton rest frame. 

In the cm- frame the 4- momenta P^, P^'', and take the form 



^b,CM 



ICM 



— {Pa.CM' 0' 0' Pa. 



CM) 



1, 0, 0, 1 



(n%/,o, 0, p3^^,) « ^(1, 0, 0,-1), 

{qo,CM, qT,CM, 0, qL,CM) , 



(39) 
(40) 
(41) 



where the simple relation between the hadron momenta and ^/s holds if the hadron masses are neglected. Note 
that without loss of generality the transverse part of the photon momentum is pointing into the a;-direction. To 
shorten the notation we will use qt = qT,CM in the following. Equations ([55 |) - ([iT|) fix the axes of the cm- frame 
according to 



qT_ 
qr 



°y,CM = ez,CM X Cx^CM , 



e-z,CM 



Pa, 



CM 



IP 



(42) 



a,CM\ 



To make the transition from the cm- frame to the CS-frame one can apply two subsequent Lorentz boosts 27|. 
In a first step one boosts along the z-axis such that the virtual photon no longer has a longitudinal momentum 
component. In a second step one boosts along the a;-axis such that also the transverse momentum of the virtual 
photon disappears. This leads to the following transformation matrix between the two frames: 



/ 



qo,cM 

-smaqo,cM 


cos a qL^CM 



-pq 

(cosa)"-^ I 





— qL,CM 

sin a qL, CM 


cosaqo^CM 



\ 



(43) 
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with 



sma = 



cos a 



(44) 



Applying the transformation matrix in (|43p to the 4-momenta cm^ ^bCM^ 1cm '-"^^ finds, in particular, that 
the hadron momenta span the xz-plane. The results are 



PafiS = {Pa,csi -sina|Pa,cs|, 0, cosa|Pa,cs|) ~ Pa,cs{^^ -sina, 0, cosa) , 

{PbfiSi -sinQ;|Pb,cs|, 0, - cosa iPfe.csl) ~ A°cs(l> -sina, 0, - cosa) 



^b,CS 



q^s = (g, 0, 0, 0) 



(45) 
(46) 
(47) 



where the energies of the hadrons in the CS-frame are given by 



pO 

^a,CS 



Pg ■ g _ \/5 / N pO Pb ■ q _ \/s / , \ 

q 2q ^ ' ' q 2q ^ ' 



(48) 



The approximate expressions in (|l5|) and (HH]) again hold if the hadron masses are neglected. Note that in this 
case one has a — a. Equations (145p . ()46p imply that the axes in the CS-frame are fixed by the hadron momenta 
according to 



1 



2 sin a V \p 



cs 



Pi 



b,CS 



a,CS\ 



\P 



b,CS\ 



ey,CS = Gz.CS X ex,CS , Sz^cS 



1 



P^ 

2 cosa \\p 



cs 



p 



b.CS 



a,CS\ 



\Pb.. 



cs\ 



(49) 



In principle there are infinitely many dilepton rest frames. Any other dilepton rest frame is related to the CS- 
frame through a 3-dimensional rotation. For instance, the frequently used Gottfried- Jackson frame (3^ . in which 
the momentum of one of the hadrons is pointing into the z-direction, is connected to the CS-frame by a rotation 
about the y-axis. 

One can readily invert the Lorentz transformation in (j43|l and find 



^ qo.CM sin a go, CM cosagi^cA/^ 

pq {cos a)^^q 

q 

yqL.CM sina qL,cM cosa9o,CJ\/ J 



(50) 



This inverse transformation is now applied to the 4- momenta of the outgoing leptons, which in the CS-frame take 
the simple form 



(l, si 



/f, _ q 



cs 



smVcs cos (pes, sm Ucs smcpcs, cost/cs 



1, - sin 9c s cos (fics, - sin 9c s sin (pes, - cos 9c s 



(51) 
(52) 



i.e., the directions of both leptons are specified by the same two angles 9cs and </)cs- This feature, of course, 
holds in any other dilepton rest frame as well. In the cm-frame the lepton momenta are given by 



/ (1 + s 



'cm - 2 



sm a sm tics cos (pes ) qo,CAi + cos a cos tics qL.CM 
qr + (cos a)~^ sin 9cs cos (pes q 
sin 9c s sin (pes q 
(l + sin a sin 9es cos (pes) qL.CM + cos a cos 9es qo.CM 



(53) 
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7'M 



( (l - sin a sin 9c s cos (pes) Qo^cai - cos a cos 9cs Il^cm \ 
qx — (cosa)"^ singes cos 0cs 9 
- sin 9c s sin 4>cs 9 
(1 - sin a sin 6'cs coshes) 9l,cm - cos a cos 9cs qo,CM j 



(54) 



By means of these momenta one can carry out the contraction of the leptonic and the hadronic tensor in the 
cm-frame. This is particularly convenient in connection with the parton model calculation in Section VI. 
We close this section with a brief discussion on the hadron spin vectors. In the cm- frame one can write 



^a,CM ^ { SaL.CM ^—-^J—^, \SaT,CM\ COS<pa,CM, \SaT,CM\ sin0a,CM, SaL,CM 



pO 



'-'b,CM 



r, \Pb,CM\ 1^ I J, lo I ■ ^ C ^b.CM 

ObL,CM ^ , \JbT,CM\ COS(pb^CM, \JbT,CM\ Sin(pb,CM: — JbL,CM 



Mb 



Mb 



(55) 
(56) 



with the longitudinal components SaL,CM, SbL,CM, and the transverse components SaT,CM, SbT.,CM- The condi- 
tion Sa = —1 implies {SaL,CM)^ + {SaT,CM)^ = 1 (and analogously for the hadron Hb). One can also write down, 
e.g., in the CS-frame in terms of longitudinal and transverse components.^ Mainly for the following reason 
we prefer, however, to work with components of the spin vectors in the cm-frame. If one has a pure transverse 
polarization in the cm- frame (in the a;2;-plane), this implies also a longitudinal polarization component in the CS- 
frame. Therefore, longitudinal and transverse polarization components can get mixed up when switching between 
both frames. Since an experimental setup and also the parton model approximation have a closer connection to 
the cm-frame than to the CS-frame it is preferable to work with cm-frame components of the hadron spin vectors. 



ANGULAR DISTRIBUTION OF THE CROSS SECTION 



By means of the general form of the hadronic tensor as derived in Section III one can now write down the full 
angular distribution of the DY cross section. Since the hadronic tensor is frame-independent this can be done, 
in principle, for any reference frame. We focus here on a dilepton rest frame because in that case the angular 
distribution takes the most compact and transparent form. Expressing the orientation of the leptons through the 
CS-angles 9cs and (j)cs (see Eqs. ([5T|) . (|52|) . and ([53| . ((54|) ) and contracting the leptonic tensor in ([5]) with the 
hadronic tensor one finds the following general form of the cross section in Eq. ([TI 



da 



d^qdn Fq^ 

' {I + COS' 9) F^j^ 

5aLf sin 26'sin0F; 



(1 - cos^ 9) F^fjj + sin 29 cos (f) F, 



uu 



sin^ 9 cos 2<j) Fi 



uu 



[Sin 
LU 



sin 20 
LU 



+ SbL 
+ \SaT 



sin 29 sin (j) F^2 '' 



-I- cos 



+ 1^, 



bT\ 



0a ( 

in 4>b(^{ 



cos 



- sin^ 9 sin 20 F\ 
SIT? 9 sin 20 F^'^l ^ 

+ (1 - COS^ 

-isin 



I F^^ + sin 29 cos F™^ + sin^ 9 cos 20 



TiCOS 20 



sin 29 sin F^'^}*" + sin^ 9 sin 20 F 



isin 20 
TU 



I Fl,rp 



I F^/rp 



+ COS( 



sin 29 sin F, 



(1 — COS' 

sin I ^. 2q^- r)^j7isin20 
UT 



sin 29 COS F^^ 



-(- sin^ 9 cos 20 F 



UT 



suT 9 sin 20 F^jj, 

+ SaL SbL ((1 + cos^ 9) fIl + (1 - cos^ 9) fIl + sin 29 cos F^°' + sin^ 9 cos 20 Fl°l ^"^^ 



* The resulting expression looks a bit more complicated because Pa,cs not pointing in the z-direction. 
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+ SaL \SbT\ COS 0b ((1 + cos^ 0) FIj, + (1 - cos^ 0) FIj, + sin 26 cos F];^' + sin^ 9 cos 20 F^^' 



COS 20 ^ 



+ sin 0b sin 26* sin Fj^jn + sin^ 9 sin 20 F^j, 



+ \SaT\S, 



bL 



cos0a((l +cos2 6')F^i + (1 - cos^ 6I)F|^ + sin 26* COS 0F^°''' + sin2 6* COS 20 F^'5^" 



[COS 2(p 



+ sin 0a ( sin 26* sin (p F -pj^ + sin^ 9 sin 2(/) i^^^^. 



■sin 20 



bT 



cos(0a + 0b) ( (1 + cos^ 9) F^j. + (1 - cos^ 0) F^j, + sin 29 cos F^^' + sin^ cos 20 F 



cos 20 
TT 



cos(0a - 0b) (^(1 + cos^ 0) F^y + (1 - cos^ 9) F|y + sin 29 cos F^'i^" ' 
sin (0a + 0b)(sin20sin0Fy^'^ + sin^ 0sin20F^'^2'^^ 
sin (0a — 0b)(sin20sin0F^'^'^ + sin^ 0sin20F^'^2'^^ | . 



sin^ 9 cos 20 F^^^ ^"^^ 



(57) 



In Eq. (|57p 48 structure functions show up which exactly matches with the number of the Vi defined in Section III. 
The structure functions again depend on the three variables Pa -q, Pb- q, and q^, i.e., Ffj^ — Fljjj{Pa • g, Pb- q, q^) 
and so on. We refrain from giving the explicit relations between the structure functions in ()57p and the Vi because 
these lengthy formulae are not needed for the following discussion. In order to shorten the notation in (1571) we left 
out indices for the angles which characterize the lepton momenta and the transverse spin vectors of the hadrons. 
There is yet another reason for omitting those indices: the form of the angular distribution in ()57p holds for 
any dilepton rest frame and not just the CS-frame. The numerical values of the structure functions of course 
change when going from one frame to another. Furthermore, note that the components of the spin vectors can be 
understood in different frames like the rest frame of one of the hadrons, the cm-frame, or a dilepton rest frame. 

In particular for the angular distribution of the unpolarized cross section different notations can be found in 
the literature (see, e.g., [37| and references therein). Here we just quote the frequently used formula 



dN 



da 



da 



1 



dn d^q dn / d^q 47r A + 3 



1 + A cos^ 9 + fisin 29 cos ( 



One readily finds 



f; 



Fi 



F2 



uu 



■sin 6* COS 20 



cos 2(p 
UU 



Fi 
^uu 



P2 
^UU 



Fi 

^uu 



P2 
^UU 



The so-called Lam- Tung relation [35|, 1 



A + 2i^ = 1 , 

which in terms of the structure functions defined in (|57p reads 

icos 20 



Puu — 2 f; 



uu 



(58) 



(59) 



(60) 



(61) 



has attracted considerable attention in the past. This relation is exact if one computes the DY process to 
0{as) in the standard coUinear perturbative QCD framework. Even at O(a^) the numerical violation of ((501 is 
small [13]. On the other hand data for N fi^ n+ X taken at CERN jSlI, [H] and at Fermilab [43| are in 
disagreement with the Lam- Tung relation. In particular, an unexpectedly large cos 20 modulation of the cross 
section was observed, and in the meantime different explanations for this phenomenon have been put forward in 
the literature 0, [iBI, H^, \^ [H, |4^, [s^l • In Ref . [sl] it was pointed out that intrinsic transverse motion of initial 
state partons might be responsible for the observed violation of the Lam- Tung relation. In the following section 
we will briefly return to this point in connection with the parton model calculation. It is also worthwhile to 
mention that more recent Fermilab data on proton-deuteron Drell-Yan do agree with the Lam- Tung relation [5l| . 

The hadronic tensor given in Section III also allows one to find the angular distribution of the cross section for 
the specific kinematical point qx — 0. Altogether, in that case one has nine independent angular dependences 



13 



and structure functions, 



da o? 



d'^qdn 



qT=o F 



em ^ 



{ ((1 + cos^ 9) + (1 - cos^ 9) F^u) + SaL SbL ((1 + cos^ 9) fIl + (1 - cos^ 9) F^ 
+ SaL \SbT\[ sin 20 cos(0 - 0b) i {Fl°^ ^ + Ft^^))+\SaT\ SbL ( sin 29 cos(0 - (ba) ^ (^^TL * + 

+ \SaT\ \SbT 



isin 0> 
TL . 



sin^ 



+ cos(0a - 0fc) ((1 + cos^ 6*) + (1 - cos^ 6*) F^y^ I . (62) 

This result was already given in Ref. using a different notation. We note that in the double-polarized sector 
the following relations hold: 

j-icos 1 j-isin I T-iCos 1 j-isin I 77icos20| j-isin2(/>| ( CK^\ 

^LT \qT.=o — ^LT lqT=0' ^TL 1,^=0 ~ ^TL lgT,=0' ^ TT lgy=o~^TT 1,^=0 ' 



All the structure functions that show up in ([57|) but not in (|62p have a kinematical zero at = 0. Notice that 
the angular distribution of the gT-integrated cross section agrees with the one in ([62]), but the corresponding 
structure functions differ numerically (see also fl3|). 

We have pointed out that our hadronic tensor provides the angular distribution of the cross section at qT — 0. In 
fact, this statement is not totally correct. The hadronic tensor in Eq. ([55]) does not generate the term proportional 
to cos(2</) — (f)a — 4'b) in ((62|) . This can be understood by taking a close look at Eq. ([3T|) : there the prefactor in 
front of h'^b 48 vanishes for qx — implying that at this particular kinematical point one is actually not allowed 
to eliminate this tensor which does generate the required cos(20 — <t>a — (t>b) term for qT = 0. One may instead 
eliminate, for instance, the tensor h'^b 47- Since the rather specific case qx — was already worked out in the 
literature we have proposed the tensor in ([55]) for the sake of symmetry. 

In Section III.D we have considered the case of identical hadrons in the initial state and the resultant constraints 
for the structure functions Vi. One can do a corresponding analysis for the structure functions defined in Eq. ([57)1 . 
The key ingredient of such an analysis is that the cross section remains the same if the hadrons are exchanged. 
Note that the exchange Ha ^ Hb also leads to the reversal of the z-direction which, in particular, implies 

^a^-^b, (l>^-4>, 9^Tr~9. (64) 

Twenty structure functions are either symmetric or antisymmetric under the exchange Pa ^ Pb- Using the 
shorthand notation of Eqs. ((5^ . (|38p one finds 





Fljuia.b), 


Fhu{b,a)^ 


FhuM, 




^-Fi;^\a,b), 


Firu''ib,a) 


= F™^2na,6), 


FhL{h,a) = 


FIlM, 


FlLib,a)^ 


FIlM, 


FriHb,a) 


= -Fri\a,b), 


Fri'Hb,a) 


^Fl°l'^{a,b), 


F^rp[b, a) — 


F^j.{a,b) , 


F^rj,[b, a) — 


F|j,(a,6) , 


F^'^{b, a) 


= -F^^'*(a,6) , 


F^°f*{b, a) 


= F™^^2^a,6), 


F^rp[b, a) — 


F^j,{a,b) , 


F^rj,[b, a) ~ 


F^rj,{a,b) , 


F^'^{b, a) 


= -F^^'*(a,6) , 


F^°f^{b, a) 


= F™^^2^a,6), 


F^^'^{b, a) 


= ~F^^'^{a,b) , 


F^%^^{b, a) 


= Ftt ) 


F^^'^{b, a) 


= F^^'*'{a,b) , 


F^'^^^ib, a) 


^-F^^^'^ia^b) 



(65) 
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The remaining structure functions fulfil the relations 



Flj^{h,a) = -F^u{a,h), F^^{h,a) = -F2^(a,6), 

a) = F™J^a, 6) , ^^^^^^6, a) = -F^°^"f'{a, b) , 

a) = ^^^'^^a, 6) , F^'^^^{b, a) = -Fj';}'^a, b) , (66) 

F^^(6,a) =F2j,(a,5), F^^{b,a) = Fl^{a,b) , 

a) = -F™^^a, 6) , F^/^^fe, a) = '^a, b) , 

F^'^l^^a) = -F2'^^a,6) , F^f'^ib,a) = Ft^'^{a,b) . 



It is of course intuitively clear that for identical hadrons relations as given in ([66]) have to exist. But one has 
to keep in mind that relative signs between the corresponding structure functions can show up. Eventually, we 
mention that ((SS)) . can also be derived from (PH)) . (|38p and the relations between the two sets of structure 
functions. 



VI. PARTON MODEL APPROXIMATION 



This section deals with the parton model description of the structure functions in Eq. ([57]) . Up to this point 
we didnot specify the external kinematics of the process. In the following we will consider the kinematical regime 
where the transverse photon momentum qt is of the order of a typical hadronic mass scale which means, in 
particular, that it is much smaller than the hard scale q. This is the region where TMDs enter the description of 
the DY process in a natural way. 

Our treatment is restricted to leading twist, i.e., to the leading order of an expansion in powers of 1/q. Mainly 
because of the potential problems of subleading twist TMD-factorization pointed out in Refs. [H, [s^ we refrain 
here from including the twist-3 case. Moreover, we neither take into account higher order hard scattering correc- 
tions nor effects associated with soft giuon radiation. For three of the structure functions such contributions were 
considered in [53| . 



A. Hadronic tensor 



The parton model description of the Drell-Yan process can be represented by the diagrams shown in Fig. [U 
where, e.g., the scattering amplitude for diagram (a) reads 



Nc ■ 9 



{Xali^t^Q) \Pa, Sa) {X,\ i'j'r''^ (0) \P,, Sb) (7^),, u{l, A) 7^ v{l\ A') 



(67) 



q c— 1 



A sum over color c and the quark flavors q is implemented explicitly in this expression. The electromagnetic 
charge of the quark, in units of the elementary charge e, is denoted by e^. A corresponding formula holds for the 
amplitude M(b) of the graph (b) . The differential cross section (fT(I|) in a dilepton rest frame is then given by 



da 



d'^qdn 8(27r) 



M, 



ib)\ 



A, A' Xa,Xb 



S'^^HPx^+Px.+q-Pa-Pb). 



(68) 



Note that there is no interference between the two diagrams in Fig. [TJ One can modify this formula by introducing 
the momenta of the active partons, ka and kb. This allows one to sum over a complete set of intermediate states 
and to rewrite the hadronic par t of the cross section in terms of fully unintegrated quark-quark correlators (see, 
e.g., Ref. [13, [13, [13, [sE, m, 113, UMl ) ■ In doing so one finds the hadronic tensor 



■^Y.^! '^"''b -^'^^ il-ka- kb) Tr[r $«(fca. Pa, SaK) $'(fcfc, Pb, Sb\nb)] + {$ ^ $} , (69) 
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where the quark-quark correlators, which depend on the full 4-niomentum of the quarks, are defined as 

/74 
e"^-^ {Pa,Sa\ i'^iO) Wdy[0, V-f (z) \Pa, Sa) , (70) 

%{h,Pb,Sb\ni,) = y^e'^-^(Pfc,^fc|V?(0)WDY[0,z|nb]V;j(z)|Pfc,5fc). (71) 

The object W denotes a gauge link operator (Wilson line) which ensures color gauge invariance of the correlators. 
We note that actually the Wilson lines cannot be derived from the diagrams in Fig. [T] They are generated, 
however, if in addition coUinear gluon exchanges between the active partons and the remnants of the incoming 
hadrons are taken into account (see, e.g., Refs. [59l. [60l. IgiI. l6^ ) . In general, the Wilson hues entering unintegrated 
parton correlators are process-dependent. For the DY process we will specify them below but already emphasize 
here their dependence on a light-cone vector Ua or Ub- Note that in Eqs. ([7n| . ([7T|) a color sum is implicit, leading 
to the factor 1/Nc — 1/3 in ([69|l . The term {<i> ^ $} in Eq. ([69|) represents the contribution of the diagram in 
Fig. [1] (b) and is obtained from the first term by interchanging the correlators. 

In the parton model initial state partons are assumed to move quasi-collinearly with respect to their parent 
hadron. Consequently, the components of the parton momenta behave like the corresponding components of the 
hadron momenta. The following estimates for the parton momenta in the DY process are valid in frames where 
the hadron Ha has a large light-cone plus-momentum and the hadron Hf, a large minus-momentum (this applies, 
in particular, to the cm- frame — see also the discussion in Section IV): 

k+^0{q), K^O{l/q), kaT^O{q^), 

k+^0{l/q), k^^O{q), hTr^O{q°), (72) 

where we use the light-cone components — {v^ ± v^)/\/2 for a generic 4- vector v. From the standpoint of 
factorization this means that $ and $ are treated as nonperturbative objects because the kinematical invariants 
ka'Pa, k^, kb-Pb, A:^ on which the correlators depend are much smaller than g^. According to ((72|) the momentum 
components k~ and k^ are small and hence can be neglected in the 5- function in Eq. (|69p . This also automatically 
implies q^ ~ fc+ and q^ ~ kj^ . The hadronic tensor then reduces to 

W'" = ^5]e2 f d%Td%TS(^\qT-LT-kbT)TY[^^<i>'^{Xa,kaT,Sa\na)j''^''{Xb,hT,Sb\^^^^ 

+ {$ ^ $} , (73) 
where we used the common DY variables 

X « ^ X,- « ^ (74) 

The transverse momentum dependent quark-quark correlators in (|73p are defined according to 

%{XaXT,Sa\na) = J ^^^^e^'--^Pa,Sa\r,{Q)mY[0,z\na]i^!{z)\Pa,Sa)\^^^,, (75) 

^%{xb,kbT,Sb\nb) = j ^^^^3 ^ e'"^-- (Pb, ^fcl (0) >Vdy[0, z\nb] ^]{z) |n, Sb)\^.^^ , (76) 

and they are obtained from the correlators in ([7n|l . ([7T|) by integrating out the respective small light-cone mo- 
mentum of the parton. We now specify the Wilson lines in the quark-quark correlators. The appropriate choice 
for the DY process is [5i,[6S[6l| 

yVDY[0, 2;|71a]|^+^Q = [0; -CO Ua] X [-00 Ua; -CO na+ Zt] X [-00 Ha + ZT ; Z^Ua + Zt], (77) 

yVDY[0,z|nfc]|^_^g = [0; ~oo Ub] X [~oonb ; -coub + zt] X [-00 Hb + zt ; z^Hb + zt], (78) 

with [a; b] denoting a straight gauge link between the positions a and b, and = (0, zt,0). The light-cone vectors 
in ([77]), jlHl) are given by 

< = -i= (1,0,0,-1), <--i=(l,0,0,l). (79) 
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Note that the diagram in Fig. [IJb) generates, e.g., the correlator <i>'(fcf,, Pf,, S'fcjn;,) which can be related to 
$'(fca, Pa, Sa\na) in Eq. (|70p by means of the parity transformation. 

We also mention that so-called light-cone divergences, which are caused by the light-like Wilson lines 
in ((771) . GHl), can be avoided if near- light-cone directions for the Wilson lines are chosen instead. For a dis- 
cussion of such divergences and other nontrivial issues concerning the precise definition of unintegrated parton 
correlation functions we refer to the recent contribution [g^I as well as references therein. 



B. Transverse momentum dependent parton distributions 

The quark-quark correlators in Eqs. ([73), ([73 can be parameterized through TMDs [H, HI, HI SI, HE Ull • A 
common and rather convenient procedure for performing such a parameterization is by specifying the traces of 
the correlators with the Dirac-matrices F = 7'^, 7^75, ia^'^'j^, 1, 275, 

$-3 [r] = i Tr [$« F] . (80) 

In the cm-frame, where the hadron Ha has a large plus-momentum, the leading (twist) traces are 
and $["^'^^5] — {1,2}), while all the other traces are suppressed in the cross section by at least one power 
of the large light-cone momentum (and consequently by one power of q). These traces then have the following 
expressions in terms of leading twist quark TMDs (see, e.g., [23. [56l|^): 

$^[7-1 = mxa,kl^) - 'TKrSiT fMi^:,^^kl^) , (81) 

$9 [7+75] = 5^^<^?J^^,fc2^) + ^_^g^^(^^,fc2^), (82) 
+ SaL ^ hiliXa, Par) + ^ ht\Xa, klr) ■ (83) 

For brevity we omitted the arguments of the correlator $. Note that the components of the nucleon spin vector 
in (l5T1) - ([55)) are understood in the cm-frame. The object represents a short form of the transverse epsilon 
tensor e , where we use the convention e '"^^ = 1. The transverse momentum dependent unpolarized quark 
distribution, helicity distribution, and transversity distribution are denoted by /i, oil, and ft-i, respectively. 
Of particular importance are also the time-reversal odd (T-odd) Sivers function [M . [65^ and Boer-Mulders 
function hj^ fss'l as they can give rise to quite interesting single spin and/or azimuthal asymmetries in hard 
semi-inclusive reactions. 

The correlator in Eq. (17^ is related to the correlator which defines, precisely in analogy to the Eqs. (|8ip- 
([83]) . antiquark distributions. For the different Dirac traces the relation reads [l^l 

[ - for7''75, 1, «75 

Since the correlator $ in (|76p is associated with the hadron Hh having a large minus-momentum in the cm-frame, 
the leading traces are now \ and $1*'^ '•'^1. Taking ([84]) into account the parameterizations can be 

directly obtained from ([5T|) 



$^[7-1 = /f(x,,fc,V) + "^l!^^^^ /lT^(^^,fcgT), (85) 



M, 



$.[7-75]= _5^^^.-J^^_42^)_:^_^g.-^(^^_fc2^)^ (86) 



5 Note that the l.h.s. in Eq. (16) of [si] should read #[»<'"^75] 
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$9 



9' h'^(x. P ) + ^hT(fc^T ■ SbT) - 



It) 



SbL ^ htlix,, kl^) - M hi\x,, kl^) . 



Mb 



(87) 



Note the respective sign change in front of the epsilon tensor which is due to the interchange of plus-momenta 
and minus-momenta. 



C. Leading spin observables 



Now we are in a position to calculate all the leading twist observables for qt q hy inserting the traces 
([ST|) -(I55 |) and ([55|) -(I57 |) into the hadronic tensor ([7S|) . We mention again that the contraction of the hadronic 
and the leptonic tensor is performed in the cm-frame where, in order to get the leptonic tensor in that frame, 
use is made of Eqs. (|54p . Since the lepton momenta contain angles in the CS-frame our final result for the 

cross section is of the form (f57|) . Here one has to keep in mind that the leading twist calculation of course merely 
provides nonzero results for part of the structure functions in (|57p . Carrying out the contraction of the tensors 
and keeping only the leading contribution in 1/q one finds 



da 



d^qdn 



^1 -I- COS' 
2 



9 



^$9 [7+] [7 ] _|_ $9 [7 + 75] (j)<3 [7 75] 

+ sin^ e ( COS 20 - 5'^5^^) + sin 24, {S'H^'^ -f t""" 

^ $} + ©(l/g) . 



75J <j)9 [vy 75J 



(88) 



To present the leading twist spin observables we will make use of the following notation for the convolution of 
TMDs in the transverse momentum space: 

C[w{kaT,hT)flh] = 



9 



j d^kaT d^kbT 5'^'^\qT - kaT - hr) w{kaT, hr) 



bT) 



(89) 



The two terms on the r.h.s. of (j89p are generated by the two diagram in Fig. [TJ For the parton model calculation 
it is convenient to introduce a number of linear combinations of various structure functions given in Eq. ()57[) : 
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if/?' 

2 \^-^TT 



cos 20 



psin 20 
^TT 



F 



cos(20— 0„+0b) 



TT 



2 V -^TT 



pCOs(20 + 0a^0fa) 1 / T-iCOS 20 y-isin 20 

^TT = 2l-^TT ^ ^TT 



F 



COs(20+0a+06) 



TT 



F 



cos 20 
TT 

cos 20 



T-isin 20 
^TT 



TT 



F 



isin 20 



TT 



Using the unit vector h 
CS-frame: 



Fuu - 



T-iCos 20 
^UU 



qT I qT one eventually finds the following leading order structure functions in the 
C [/i A] , (91) 

C -^T){h-kbT)-KT-kbT ^ J ^g2) 
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(93) 
(94) 
(95) 
(96) 
(97) 
(98) 
(99) 
(100) 
(101) 
(102) 
(103) 
(104) 
(105) 
(106) 
(107) 
(108) 
(109) 
(110) 
(111) 
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(112) 
(113) 



f: 



TT 



'4:{h ■ kaT) {h ■ Ut) [2(K • Kt) {h ■ Ut) - kaT ■ k 



bT\ 



+ 



'^aTHT 



2kl^{h-kt,Tf -2kl^{h-kaT) 



(114) 



We close this section with a number of comments. 



• The structure functions depend on the variables {xa,Xb,qT)- Instead of using qt one may also work with 
the transverse momentum of one of the hadrons in the CS-frame. 

• One finds nonzero contributions for 24 out of the 48 structure functions defined in Eq. (j57p . This also means 
that exactly half of the structure functions are of subleading twist for the kinematical region qx q we are 
interested in here. 

• The leading twist parton model calculation containing T-even effects was first carried out in Ref. fs^J , while 
T-odd effects were investigated in iSS] . We obtain the same number of nonzero structures identified in those 
articles, though we do not agree with certain angular dependences given in (s^ . 



• Our results are for the structure functions in Eq. ([57]) with the lepton angles understood in the CS-frame, 
and the components of the hadron spin vectors in the cm-frame. Note that the expressions would be exactly 
the same for structure functions defined in the Gottfried-Jackson frame, because differences between those 
two dilepton rest frames are only of O^qx/q)- 

• For identical hadrons in the initial state the results in Eqs. (PT|) - (|114p satisfy the model-independent con- 
straints listed in (|65p and (|66p . In particular, we point out that the parton model result 



TT 



{Xbj Xa) — Fr 



COs(2lp — lpa+<pb) 

TT 



(115) 



has a model-independent status. It is worthwhile to mention that, by means of charge conjugation, in 
the case of proton- antiproton DY one also finds symmetries for structure functions (like Firij{xa, Xb, qT) = 
Fuu{xb, Xa, qT)), and relations between various structure functions. In particular, when studying single spin 
effects one can obtain the same information by either polarizing the proton or the antiproton. 



• If the cross section is integrated upon qx only three structure functions {F^jj, F^j^, Fr 
Neglecting hadron masses one obtains 

da s da 
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dxa dxb d^ 2 dq~^ dq 
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+ \SaT\ \SbT\ sin^ e cos{2^ -^a- <l>b) ^ 6^ (/if (a;,) hfixb) + hf{xa) hlixb)) } . (116) 



Further integration upon the solid angle Q provides 



da 



dxn dxb 



9q2 



{fiMf!ixb) + f!ixa)f!ixb) 
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(117) 
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Note that the term containing the transversity dropped out. 

• For the gT-dependent cross section ah chiral-odd parton distributions disappear after integrating out the 
azimuthal angle (p. On the other hand, all the chiral-even effects survive this integration. 

• The large number of independent structure functions — for instance 16 for identical hadrons in the initial 
state — indicates the high potential of the polarized DY process for studying TMDs. Therefore, this process 
has also a certain advantage over semi-inclusive DIS (if in that reaction polarization of the initial state Icpton 
and hadron are exploited) where eight leading twist structure functions exist [23l. |2^. being just sufficient 
to map out, in principle, all the eight leading twist TMDs. 

• As already pointed out in Section V data on tt" ^ /i+ X [U, IH, [i^ show a rather large cos 2(j) 
dependence of the unpolarized cross section which cannot be explained by collinear perturbative QCD. 
However, if intrinsic transverse parton motion in the initial state is taken into account the Boer-Mulders 
function h-^ contributes to the cos 2(f) term according to (|92p which may explain the observed violation of 
the Lam- T ung relatio n Isp. This findin g st imulated a lot of phenomenological work on this subject [HI, [H, 

• Of particular interest is also the transverse single spin effect given by F^jj in Eq. jMI) or FIjj. in (|98l). 
Both structure functions contain the Sivers parton distribution which was predicted to have the opposite 
sign in DY as compared to semi- inclusive DIS [H, [8l|, [s^l- As the sign reversal is at the core of our 
present understanding of transverse single spin asymmetries in hard scattering processes an experimental 
check of this prediction is of utmost importance. Theoretical work on the Sivers effect in DY can be found 

in pSHIslIsi,!!^,^!!^. 

• The expected sign reversal of T-odd TMDs can also be investigated through the structure func- 
tions Fj^ji^ '^"^ in (|96p or F^'J^^"^ in ([U^ in which the Boer-Mulders function enters (sec also 
Refs. [illMI^)- 

• A phenomenological study of the structure functions in (|109|) . piOp was carried out in [90l |. 

VII. SUMMARY 

We have presented a formalism for dilepton production from the collision of two polarized spin-^ particles. To 
this end we have derived in a first step a general expression for the hadronic DY tensor. This tensor consists 
of 48 basis elements, and each basis tensor is multiplied by a scalar function (structure function). In order to 
ensure electromagnetic gauge invariance of the hadronic tensor we have made use of an elegant projection method 
proposed in [1^. In general, our treatment completes earlier work [l3,[3l|. The double polarized case, which is 
the most challenging part, was studied before only for the specific kinematical case = [13] ■ 

The result for the hadronic tensor allows one to obtain the general angular distribution of the cross section for 
any reference frame. In this work we have focussed on a dilepton rest frame where the angular distribution takes 
the most compact form and shows a high degree of symmetry. We repeat here that the angular distribution as 
given in Eq. (|57p . which represents a central result of our work, holds for any dilepton rest frame. 

Our analysis is supplemented by a parton model calculation of the polarized DY reaction (see also [H, Issj). 
For this part of the work we concentrated on the kinematical situation where the transverse momentum of the 
dilepton pair is much smaller than its invariant mass. This region is the realm of TMDs which are currently under 
intense investigation both from the experimental and the theoretical side. 

We reemphasize that the polarized DY process has a high potential for studying TMDs which contain important 
information on the nonperturbative structure of hadrons. Moreover, polarized dilepton measurements can provide 
us with a crucial and highly nontrivial check of QCD-factorization [59| . In addition, one can systematically study 
different resummation techniques [oij . [92 . 93 ] in an unprecedented way. Consequently, there is sufficient reason 
for looking forward to the first polarized DY data. 
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